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The response (e.g., dose-rate) of an isotropic detector to primary radiation from a 
finite plane source may, for points less than a mean-free-path distant, be evaluated as the 
sum over an infinite series, This series is derived by expanding the exponential dependence 
exp(— Mor) as a power series in nor and integrating, numerically or analytically, each term over 
the source array. If the medium is effectively infinite and homogeneous, the scattering 
properties of this medium can be characterized by a point isotropic source buildup factor. 
Buildup factor data are often approximated by formulas having simple analytic dependences 
on iJLor and numerical parameters independent of mo^. Any such set of parameters can be us(kI 
to generate a set of weight functions bn for an infinite scries buildup factor representation 

exp i + fior) X) &«• (-Mor)«M! 

71 = 

which can be used with the above primary radiation series solution to give, instead, a 
series solution including both primary and scattered radiation. Tables of bnS (air-dose 
buildup only) derived from coefficients of cubic polynomials fitted to the Goldstein- Wilkins 
data (NYO-3075) are given for <n < 13 at primary photon energies of 0.5, 0.7, 1.0, 1.5, 2.0, 
3.5, 5.5, 7.5, and 9.5 Mev in water, Al, Fe, Sn, W, Pb, and U. The; method is applied to give: 
(1) a comparison with dose-rates measured by a detector separated from a Co-OO rectangular 
plaque food irradiator by a layer of steel and a layer of water and (2) dose-rate profiles at 
constant heights across a cleared circular area in a Co-60 infinite plane source in air. 



1. Introduction 

A number of methods for evaluating radiation 
fields from distributed sources have been developed 
as part of a general program at the National Bureau 
of Standards [1, 2]^ and otlier agencies. Some of 
these methods, which provide easy and rapid solu- 
tions to radiation shielding problems, rely on the 
approximate proportionality between dose-rate and 
the solid angle subtended by a uniformly radiating 
surface. This approximation has been shown to be 
quite good if the material interposed between source 
and detector is of the order of a mean free path in 
thickness [3]. 

For situations involving much less than a mean 
free path of interposed material, as in the case of a 
detector separated from a surface deposition of 
radioactive material by only a thin roof or a few feet 
of air [4, 5], the power -series method developed by 
Sievert [6] for circular disk sources has been extended 
[7, 8] to include rectangular sources and to take into 
account point-source buildup factor data in poly- 
nomial form [9]. 



nVork sui)i)()rl(Ml jointly hy Uw OlTnv, of Civil Defeuse (DOD), the Navy 
Bureau of Yards and Docks, and tiic DclViisc Aloiiiic Supi»()rt Agency. 

1 Figures in brackets indicate tiic literature references at the end of this paper. 



In addition, it has been reported [lOJ that buildup 
data from a number of other analytic foi'inulations 
[11-15] can also readily be transformed for direct 
use in a Sie vert-type power-series solution. The 
present work constitutes a published account of the 
material reported in [10] and includes a quantitative 
comparison of the buildup data foj'inulations con- 
sidered [9, 11-15] for Co-C)0 gamma rays in aluminum 
for zero to ten mean-free-patii penetrations, a tabic 
of buildup data, ti'aiisformations, a sample table of 
transformed buildup data, and sample calculations 
and results. 

Besides fallout radiation studies [5], the power 
series method is applicable, or has already been ap- 
plied [16] to such diverse distril)uted source studies 
as: 

(a) Tracer studies involving an elongated rec- 
tangle as the initial distribution of radioisotope- 
tagged sand in a beach-erosion sand-drift experi- 
ment [16], 

(b) food irradiator design studies involving plaque 
sources [8, 17, 18], and 

(c) calculation of the dose distribution in tissue 
adjacent to disk or other shaped radioactive appli- 
cators for treatment of malignancies [6, 19]. Also, 
buildup data in the form of the series coefficients 
presented here are directly applicable to a method 
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of barrier shape optimization described by Jain and 
Sharma [20] who use, in effect, the first three terms 
of this series for their analysis. 

2. Buildup Factor: Definition and Discussion 

As in references [9] through [15], the source and 
detector environment is ideahzed throughout this 
paper as a homogeneous infinite medium, shown 
schematically in figure 1. In such an idealized 
situation the scattered radiation reaching the de- 
tector can be taken into account by means of a 
point isotropic (PTI) source buildup factor [21] 



£pTi(So, Z,r) = (Z}»+Z>')/Z)» 



(1) 



which depends on the effective atomic number Z 
characterizing the medium and the distance r be- 
tween the detector and a monoenergetic point 
source of photons of energy ^ E^. 

The quantities D^ and ^Mn (1) refer to the un- 
scattered and the scattered components, respec- 
tively, of the radiation received by the detector. 
For a given physical situation, D^ and D^ (and 
hence Bp^iiEo, Z, r)) can assume a variety of nu- 
merical values depending on what kind of detector 
response (e.g., photon number flux [22], energy 
current [22] at an interface, energy dissipation in a 
medium- or air-equivalent detector, etc.) is used to 
describe the radiation field. The buildup factor 
data in this work correspond to a detector which 
measures radiation ^^/:x in terms of energy dissipation 
in air. In this case Bpti{Eo, Z, r), D^ and D^ in (1) 
are the same as the quantities Ba, h^^^ and la"^^ 
defined in [21]. 



2 Here, and in what follows, the source is considered to be either monoenergetic 
or an effectively monoenergetic component of a photon energy distribution over 
which the final results must be summed. 




Figure 1. Point isotropic source geometry. 



The total detector response Z) to a point isotropic 
source is 

D=D'+D^=D'B^r,,(Eo,Z,r) 

= (/:/47r)[exp (-iJior)/r']B^Ti(Eo, Z, r) (2) 

in which fjLo[= jjl^Eq, Z)] is the narrow-beam attenua- 
tion coefficient [23] at photon energy Eq for the 
medium. The source strength factor k, for the de- 
tector response-type {D=Id) specified above is 

k=nXEQXiJien(Eo, Z) X (dimensional factors) (3) 

where n is the total number of photons of energy Eq 
emitted per second by the source, fXen(Eo, Z) is the 
energy absorption coefficient [24] for air,^ and the 
^^dimensional factors'' provide the units [e.g., (rads/ 
hr)/curie] desired in the final answer. 

3. Analytic Formulations for ^ptiC^/o, Z, r) 

For applying buildup data to distributed source 
problems, the dependence of J5pTi(£'o, Z,r) on r is 
usually approximated by some analytic function. 
This not only provides some economy in data tabu- 
lation, but the function can be so chosen that its 
presence entails little or no additional complication 
in the integration 



Z>=(<T/47r) £[exp(-Mor)/r^]SpTi(£o, Z, r)dS 



(4) 



of (2) over a distributed source S with differential 
elements dS, as schematized in figure 2. The con- 
stant factor in (4) 

CT = dk/dS (5) 



3 Note, however, that example 1 in section 7 requires an answer in terms of 
energy dissipation in ivater {D=^It>) rather than in air (D=Id) , so, that m en (EofJiiO) 
is used in the computation in 7.1, rather than Men (Eo, air) corresponding to the 
data in table 3 in section 5. 




Figure 2. Finite plane isotropic source geometry. 
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is the source strength per unit of source area or 
volume, with /: defined as in (3). 

The following functions (a) to (f) have been used 
to represent point isotropic source BpTi{f^Jo, Z, r) data 
in order to obtain analytic solutions to (4). 

(a) Linear approximation [11 J (one parameter): 



where 



BpTiiEo, Z,r)==l +a^or 



«=(Mo— Men)/Alej 



(6) 



(7) 



(b) Exponential-linear empirical fit [12] (two 
parameters): 



BpTi(Eo, Z, r)=l+aiJior exp (— &mo^). 



(8) 



(c) Exponential-polynomial fit [13] (one or three 
parameters): 

5pTi(^o, Z, r)=cxp (+nor)(br+l) (9a) 

Of, for larger thicknesses: 

5pTi(So, Z, r)=exp i+^or)(ar'+br+c). (9b) 

(d) Low order polynomial fit [9, 25] (three 
parameters) : 

i?PTi(So, Z, r) =l+^iMa^+&(/xor)2+/33(Mor)^ (10) 

(e) Sum-of-exponentials fit [14] (three parameters) : 
BpTi(Eo,Z^7')=Aexp {—aiHo^) + (l—A) exp (—aoMr) 

(11) 



(f) Exponential-linear moments calculation result 
[15] (four parameters) : 

BpTiiEo, Z, r)-l+Mo^{^i2?i exp \{\-B,)ix,t\ 

+^2fi2e.xp[(l-i?2)Mor]}. (12) 

Table 1 contains data calculated at /xor=l, 2, 4, 7, 
and 10 mean free paths according to approximation 
(a) above, empirical fits (b) to (e) and direct-result 
parameters (f), in addition to earlier cak'ulated 
results from NYO-3075, [26] (g), all for ^'exposure 
dose'' (air-equivalent detector) in aluminum and 
interpolated to 1.25 Mev. 

The column headed e/ beside each column of 
calculated data contains percent deviations from the 
recent Berger-Spencer [15] results (f). 



^ ^ B-B (Berger-Spencer) ^^^QQ^ 
^ B (Berger-Spencer) 



(13) 



and the column headed e^ contains percent deviations 
from the Goldstein-Wilkins results (g) to which (b) , 
(d), and (e) were fitted, 



B—B (Goldstein-Wilkins) ^ 

^'~ /^(Goldstein-Wilkins) ^^^^ ^^' 



(14) 



As seen in table 1, the Goldstein [11] linear 
approximation (a), above, involves neither fitting 
nor transport theory calculations. It is very useful 
for penetrations up to one or two mean free paths, 
and gives results within a factor of two up to ten 
nu^an free paths in aluminum. 



Table 1. The point isotropic source buildup factor {"exposure dose": air-equivalent detector), Bvti{Eq,Z), for 1.25-Mev gamma rays 
in aluminum evaluated from five different analytic approximations (a to e) at 0, 1, 2, 4, 7", and 10 mfp, compared with the 
corresponding moments-calculated analytic (/) and numerical (g) results 





(a) Goldstein [11] 1-parameter 
linear approximation (6) 


(b) Chilton et al. [12] 2-pa- 

rameter exponential-linear fit 

(8) 


(c) Leshchinskii [13] l-parani- 
eter exponential-linear fit 
(9a) 


nor 


J5PTI 


«/* 


^.** 


-BPTI 


</ 


^g 


BPTI 


«/ 


^g 




1.00 
2.06 
3.01 
5.04 
8.04 
11.0 


% 

7.9 
-1.0 
-13 
-27 
-36 


% 

7.3 
-2.0 
-13 
-29 
-37 


1.00 
1.97 
3.06 
5.69 
11.0 
18.3 


% 



3.1 

0.7 

-2.2 



6.4 


% 


2.6 
-0.3 
-2.2 
-3.5 

4.6 


1.00 

1.80 

2.41 

-18.9 

-1490 

-52100 


% 


-5.8 
-21 





1 
2 


\ 

-6.3 
-21 


4 






7 






10 













(e) Taylor [14] 3-parameter 


(f) Berger-Spencer [15] 4-pa- 


(g) Goldstein-Wilkins [26] 




(d) Capo [9] 3-parameter poly- 


sum-of-exponnetials fit (11); 


rameter exponential-linear 


moments calculation re- 




nomial (cubic) fit (10) 




Lakey [28] data 


moments calculation re- 
sults (12) 


sults; Al data in tabulated 
form only 


Mor 


^PTI 


«/ 


^g 


BpTl 


«/ 


^g 


J5PTI 


^PTI 






% 




% 




% 


% 









0.997 


-0.3 




-0.3 


1.00 








1.00 


1.00 


1 


1.93 


1.0 




0.5 


2.10 


9.0 


9.4 


1.91 


1.92 


2 


3.05 


0.3 




-0.7 


3.26 


7.2 


6.2 


3.04 


3.07 


4 


5.81 


-0.2 




-0.2 


5.76 


-1.0 


-1.0 


5.82 


5.82 


7 


11.1 


0.9 




-2.6 


10.2 


-7.3 


-11 


11.0 


11.4 


10 


17.6 


2.3 




0.6 


15.7 


-8.7 


-10 


17.2 


17.5 



={[BpTi--BpTi(Berger-Spenccr)]/BpTi(Berger-Spcncer) j X100%. 
={[5pTi-BpTi(Goidstein-Wilkins)]/i^PTi(Goldstein-Wilkins)}X100%. 



293 



The two-parameter exponential-linear empirical 
fit (b) of Cliilton, Holoviak, and Donovan [12] is 
surprisingly good for penetrations up to ten mean 
free paths. This formulation, similar in form to 
the calculated result (f) of Berger and Spencer 
[15], can thus provide economical storage and easy 
application of buildup data over a wider range of 
penetration thicknesses. Reference [12] contains 
parametric data for aluminum only, which was the 
deciding factor in selecting this as the reference 
substance for table 1. 

Leshchinskii [13] supplied data, for the reference 
situation of Co-60 gammas in aluminum, only for 
the one-parameter exponential-linear fit (9a), in- 
tended for penetrations of to 10 cm (0 to 1.49 
mean free paths), although for Co-60 gammas in 
water he did use the three-parameter fit (9b) for 
penetrations of 2 to 70 cm (0.127 to 4.44 mean 
free paths). In each case the fit is within about 5 
percent up to the specified upper limit and then 
departs wildly. However, this formulation does 
have the advantage that the exponential dependence 
exp (-/do'T) is canceled out of the distributed source 
integral (4), and it will be seen that the power-series 
solution of (4) [8, 10] is, in effect, a generalization 
of the Leshchinskii result. 

The slight departure of the Capo data (d) from 
the condition Bfti(Eo, Z, r)=l at jjlqT^O results 
from additional cross fitting as a polynomial in 
energy or reciprocal energy, in which it is difficult 
to preserve this normalization. Tliis discrepancy 
is not large enough to detract from the usefulness 
of the data, however, and in this specific comparison 
(1.25 Mev photons in Al) the departures e/ and €g 
from calculated data (f) and (g) are seen to be 
less, over the given range of juoV values, than cor- 
responding departures of (a), (b), (c), and (e). 

In addition to the empirically fitted Capo tabula- 
tion, Pi polynomial coefficients for Co-60 and 
Cs-137 gamma rays in concrete and water were 
directly calculated by Spencer [1] using the moments 
method and were punched on IBM cards for use 
in further calculations. In [1] the Aj coefficients 
in (Al), page 80, are identical with 13 i coefficients 
if (a) Aj refers to the zeroth Legendre harmonic 
(/=:0), and (b) Aj refers to integral dose data rather 
than to the differential energy data indicated by 
the dependence on Ei in (Al). 

The sum-of -exponentials three-parameter fit (4) 
of Taylor [14] is perhaps the most well-known and 
widely used buildup factor formulation in current 
use [27]. The only published parametric data for 
aluminum seem to be those of Lakey [28], and these 
were used for the comparison in table 1. 

The exponential-linear four-parameter formula (f) 
of Berger and Spencer [15] represents direct calcu- 
lated results by the moment method [21]. Parametric 
data in [15] are tabulated for 0.0341- to 10.22-Mev 
gamma-rays in aluminum and concrete, and addi- 
tional formulas are supplied for applying these build- 
up parameters to such geometries as the isotropic 
disk source (on-axis) and the isotropic spherical 



surface source. The buildup data in table 1, obtained 
by evaluating (12) using parameters interpolated to 
1.25 Mev, agree within 4 percent with the earlier 
moments calculation [21] results of Goldstein and 
Wilkins [26], (g), to which the formulas (b), (d), and 
(e) were fitted. 

4. Power-Series Solution for Unscattered 
Radiation From a Distributed Source 

It has been shown [6, 7] that for finite plane sources 
of simple shape, a solution of (4) for D^ (i.e., omitting 
the factor 5pTi {Eq, Z, r) within the integrand) maybe 
obtained in the form of an infinite series 



D'-- 



- ((7/47r) X) g.n (geom) • (/xqx)^ 

n=0 



(15) 



which converges quickly when the detector is less 
than a mean free path from the source. The dimen- 
sionless ^^(geom) coefficients in (15) depend only on 
the shape of the detector-source geometry, and are 
derived from a power series expansion of the 
exponential in (4) 



exp {-nor)=J2 (-^iorr/nl (16) 

71=0 



such that 



qn (geom)= f {-r/x)''dS/{r' • n!). (17) 



The quantity ijloX, also dimensionless, is the mean- 
free-path thickness of material between source and 
detector. 

Formulas and tables of ^,, (geom) coefficients (tables 
include values for 0<7^<9) have been generated for 
a rectangular source [7] and for a circular disk source 
with detector off-axis [8]. As a check, it was noted 
[8] that for a detector on-axis over a circular disk 
source the series (15) reduces to the familiar exponen- 
tial-integral solution, in which the exponential 
integrals Ei (x) are replaced by power series in x. 

5. Adaptation of Analytic B^ti {Eo, Z, r) For- 
mulations to the Distributed Source Power- 
Series Solution 

In the series solution (15) for the unscattered 
component D^ in section 4, the material between 
source and detector can either be distributed through- 
out the intervening space, or can be arranged as one 
or more uniform slabs or sheets lying parallel to the 
source-plane. In this section, on the other hand, 
the source and detector are assumed embedded, as 
discussed in section 2, in an infinite homogeneous 
medium. 
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ReplaciMs>: ilio polynomial in the Leslichinskii for- 
mula (91)) l)y a ])ower series we obtain the buildup 
facto I' fo I'll 111 la lion 



SpTi(£o, Z, /O=oxp (+Mo/') S hn{E,, Z) . (-Mor)>! 

7i = 

(18) 

reported in [10], from which the coefScients* hn {Eq, Z) 
can be used with the series solution (15). The only 
change in (15) is the multiplication of each series 
term g„(geoni)-(/xoCc)" by the additional factor 
bn(Eo, Z) to give the total dose-rate solution 



D=D'+D^^{<rl4.ir) S K{E,, Z) . g,(geom) . {^,xY. 

71 = 

(19) 

Equations (18) and (19) are, of course, of little 
practical value if an infinite number of bn{Eo, Z) co- 
efficients must be fitted and tabulated. However, by 
expanding as power-series all exponentials appearing 
in the analytic formulations (a) to (f) in section 3, 
then equating coefficients of like powers of fx^r, a 
transformation for generating bn{E^, Z) coefficients 
from the parameters of any of these formulations may 
be obtained. Table 2 presents these ''parasite'^ trans- 
formations, and includes (LS) for comparison. 



4 The simpler notation bn(Eo,Z), used throughout this work, is ecjuivalent to 
d'n {Eo, Z) used in [8] and [10]. 



Tlie surprisingly simple transformation (e) of 
Taylor [14] and Lakey [28] parametric data results 
because inclusion of a Bpti{Eo, Z, r) of this form into 
(4) produces no change in the analytic form of the 
integraiul. It can also be seen in table 2 that the 
Goldstein [11] linear approximation (a) is a special 
case of the polynomial fit (cl) used by Capo [9], and 
that the four-parameter calculated result (f) of 
Berger and Spencer [15] is identical in analytic form 
with the two-parameter empirical fit (b) of Chilton 
et al. [12]. The similarities are reflected in the 
respective transformations in the last column. 

To illustrate the use of table 2, transfonnation (d) 
was used to generate the ^^.(^o, Z) values for <7i < 13 
in table 3 from ''dose'' (air-equivalent detector) 
buildup data for water, Al, Fe, Sn, Pb, W, and U in the 
form of ^i{EQ, Z) polynomial coefficients. The input 
fiiiE^, Zys were taken from reference [9], tables IB, 
4B, 7B, IOC, 13B, 16C, and 18B wliich in turn were 
fitted by least squares to data calculated by Gold- 
stein and Wilkins [26] using the moments method. 
Of the 15 gamma-ray energies from 0.4 to 9.5 Mev 
presented in [9] for reactor ''operating'' and "shut- 
down" shielding calculations, the 10 energies £'o=0.5, 
0.7, 1.0, 1.5, 2.0, 2.5, 3.5, 5.5, 7.5, and 9.5 Mev were 
arbitrarily selected. 

'J'lie appearance in table 3 of b^iEo, Z) [=^q{Eo, Z) 
= B(Eo, Z,r)iiit iJLor=0} values which are not unity 
follows from the discussion of the Capo data in 
section 3. 



Table 2. — Summary of analytic formulations (a to f) for the point isotropic source buildup factor Bpti(/^o, Z) {column 4) 

The power-series fornuilation, in which the " »" for number of parameters (column 3) may be replaced by 1, 2, 3, 4, or N, depending on which of the formulations 
(a to f) is used for generating the bniEo, Z) coellicients, has l3een added for comparison. The last column contains the "parasite" transformations, of which (d) 
was used to generate the sample bniEo, Z)'s in table 3 from the Capo [9] tables of fii polynomial coefficients 



Formu- 
lation 


Reference 


Number of 
parameters 


BpTiiEo, Z, r), point isotropic source 


EqNo. 


bniEo, Z); 71=0, 1, 2, . . . 


(a) 


Goldstein [11] 




1 


l+«Mo?- 


(6) 


l—na 


(b) 


rChilton,Holoviak, 
IDonovan [12] 




2 


l+a/i/ exp i-bfi^r) 


(8) 


1-a S i C]) &-1 
i = l ^^^ 


(e) 


Leshchinskii [13] 




1 
3 


exp(+iUor)(6r+l) 
exp i+n^r) (ar2+6r+c) 


(9a) 
(9b) 


bo=l; bi = -bJfiQ; bn>2=0 

bo=c; bi = -b/n-, b2=2a/^o2; 
" bn>z=0 


(d) 


Capo [9] 




3oriV 


N 

S fiiiEo, Z)(n,Ty 
i=0 


(10) 


n<N 
2 i-mnmn-iy.]fiiiEo, Z) 
1 = 


(c) 


/Taylor, [14] 
iLakey [28] 




3 


A exp (-ai/xor)+(l-^) exp (-a2fiQr) 


(n) 


A(l+ai)''+(l-yl)(l+a2)" 


(0 


Berger, Spencer [15] 




4 


l-\-HQr[AiBi exp [-{Bi-Dfi^r] 
+yl2B2exp[-(B2-lW]} 


(12) 


l-.S i(^>) [A^B^iB,-ly-l 
■^AiB2iBi-\y-^ 


Tower 

series 


This work 




00 


exp(+Mo'') 2 bniEo, Z)'(-fis)nln\ 
n = 


(18) 


bniEo, Z) 
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Table 3. "^Coefficients bn(Eo; Z) of the power-series buildup factor formulation {18) derived from the Capo [9] ^i polynomial coeffi- 
cients according to transformation (d), table 2. 

These 6„ (Eo, Z) coefficients may be used directly in the distributed source series solution (19) with geometry coefficeints Qn (geom) (17) such as are given in references 
[7] and [8] for rectangular or ofF-axis disk sources. 

6n(£o, H2O) 



\ So(Mev) 












n \^ 


0.5 


0.7 


1.0 


1.5 


2.0 





1.0011 (0) 


3.0008 (0) 


9.9678 (-1) 


9.9480 (-1) 


9.9542 (-1) 


1 


6.4656 (-2) 


5.0643 (-2) 


9.7950 (-3) 


6.6752 (-2) 


1.5973 (-1) 


2 


2.3708 (-1) 


-2.4382 (-1) 


-6.3441 (-1) 


-7.1188 (-1) 


-5.9935 (-1) 


3 


1.4491 (0) 


1.0713 (-1) 


-9.2998 (-1) 


-1.3348 (0) 


-1.2776 (0) 


4 


3.G313 (0) 


1.0932 (0) 


-8.7105 (-1) 


-1.7959 (0) 


-1.8709 (0) 


5 


6.7143 (0) 


2.7040 (0) 


-4.5177 (-1) 


-2.0887 (0) 


-2.3751 (0) 


6 


1.0629 (1) 


4.9293 (0) 


1.3372 (-1) 


-2.2071 (0) 


-2.7859 (0) 


7 


1.5306 (1) 


7.7588 (0) 


3.4913 (0) 


-2.1448 (0) 


-3.0992 (0) 


8 


2.0675 (1) 


1.1182 (1) 


3.0268 (0) 


-1.8955 (0) 


-3.3109 (0) 


9 


2.6668 (1) 


1.5189 (1) 


4.9461 (0) 


-1.4529 (0) 


-3.4167 (0) 


10 


3.3216 (1) 


1.9769 (1) 


7.2551 (0) 


-8.1092 (-1) 


-3.4124 (0) 


11 


4.0248 (1) 


2.4912 (1) 


9.9595 (0) 


3.6838 (-2) 


-3.2940 (0) 


12 


4.7695 (1) 


3.0608 (1) 


1.3065 (1) 


1.0966 (0) 


-3.0571 (1) 


13 


5.5489 (1) 


3.6846 (1) 


1.6578 (1) 


2.3746 (0) 


-2.6977 (0) 


\j5;o(Mev) 












n \. 


2.5 


3.5 


5.5 


7.5 


9.5 





9.9663 (-1) 


9.9899 (-1) 


1.0022 (0) 


1.0042 (0) 


1.0054 (0) 


1 


2.4502 (-1) 


3.7542 (-1) 


5.2925 (-1) 


6.1392 (-1) 


6.6696 (-1) 


2 


-4.6449 (-1) 


-2.3628 (-1) 


4.9890 (-2) 


2.1199 (-1) 


3.1471 (-1) 


3 


-1.1293 (0) 


-8.3515 (-1) 


-4.3613 (-1) 


-2.0236 (-1) 


-5.2252 (-2) 


4 


-1.7467 (0) 


-1.4203 (0) 


-9.2912 (-1) 


-6.2983 (-1) 


-4.3481 (-1) 


5 


-2.3140 (0) 


-1.9906 (0) 


-1.4294 (0) 


-1.0711 (0) 


-8.3385 (-1) 


6 


-2.8287 (0) 


-2.5453 (0) 


-1.9372 (0) 


-1.5270 (0) 


-1.2503 (0) 


7 


-3.2881 (0) 


-3.0834 (0) 


-2.4529 (0) 


-1.9981 (0) 


-1.6849 (0) 


8 


-3.6895 (0) 


-3.6040 (0) 


-2.97^7 (0) 


-2.4851 (0) 


-2.1387 (0) 


9 


-4.0302 (0) 


-4.1060 (0) 


-3.5090 (0) 


-2.9888 (0) 


-2.6124 (0) 


10 


-4.3077 (0) 


-4.5885 (0) 


-4.0501 (0) 


-3.5099 (0) 


-3.1071 (0) 


11 


-4.5193 (0) 


-5.0506 (0) 


-4.5980 (0) 


-4.0491 (0) 


-3.6235 (0) 


12 


-4.6623 (0) 


-5.4914 (0) 


-5.1597 (0) 


-4.6071 (0) 


-4.1626 (0) 


13 


-4.7340 (0) 


-5.9099 (0) 


-5.7288 (0) 


-5.1846 (0) 


-4.7253 (0) 







bniEo, 


Al) Z=13 






\Eo(Mev) 












\ 


0.5 


0.7 


1.0 


1.5 


2.0 


n \ 















1.0001 (0) 


9.9999 (-1) 


9.9819 (-1) 


9.9641 (-1) 


9.9640 (-1) 


1 


-1.3070 (-1) 


-2.5215 (-1) 


9.0059 (-2) 


2.1928 (-1) 


2.7129 (-1) 


2 


-7.8756 (-1) 


-1.2030 (0) 


-5.6633 (-1) 


-4.0565 (-1) 


-3.6339 (-1) 


3 


-9.8747 (-1) 


-1.8585 (0) 


-9.6259 (-1) 


-8.6972 (-1) 


-9.0190 (-1) 


4 


-7.4740 (-1) 


-2.2242 (0) 


-1.0904 (0) 


-1.1643 (0) 


-1.3385 (0) 


5 


-8.4323 (-2) 


-2.3060 (0) 


-9.4128 (-1) 


-1.2806 (0) 


-1.6674 (0) 


6 


1.1545 (0) 


-2.1097 (0) 


-5.0695 (-1) 


-1.2101 (0) 


-1.8828 (0) 


7 


2.4430 (0) 


-1.6409 (0) 


2.2099 (-1) 


-9.4414 (-1) 


-1.9791 (0) 


8 


4.2732 (0) 


-9.0553 (-1) 


1.2509 (0) 


-4.7400 (-1) 


-1.9505 (0) 


9 


6.4586 (0) 


9.0729 (-2) 


2.5912 (0) 


2.0895 (-1) 


-1.7912 (0) 


10 


8.9821 (0) 


1.3421 (0) 


4.2502 (0) 


1.1134 (0) 


-1.4954 (0) 


11 


1.1827 (1) 


2.8321 (0) 


6.2364 (0) 


2.2479 (0) 


-1.0575 (0) 


12 


1.4976 (1) 


4.5870 (0) 


8.5580 (0) 


3.6213 (0) 


-4.7163 (-1) 


13 


1.8412 (1) 


6.5691 (0) 


1.1224 (1) 


5.2421 (0) 


2.6792 (-1) 


\^c(Mev) 












n \. 


2.5 


3.5 


5.5 


7.5 


9.5 





9.9700 (-1) 


9.9849 (-1) 


1.0008 (0) 


1.0022 (0) 


1.0032 (0) 


1 


3.2250 (-1) 


4.2042 (-1) 


5.6586 (-1) 


6.5832 (-1) 


7.2035 (-1) 


2 


-2.9618 (-1) 


-1.3354 (-1) 


1.3872 (-1) 


3.7355 (-1) 


4.4436 (-1) 


3 


-8.5539 (-1) 


-6.6191 (-1) 


-2.8031 (-1) 


-1.1564 (-2) 


1.7550 (-1) 


4 


-1.3515 (0) 


-1.1632 (0) 


-6.9091 (-1) 


-3.3714 (-1) 


-8.5906 (-2) 


6 


-1.7808 (0) 


-1.6359 (0) 


-1.0928 (0) 


-6.5622 (-1) 


-3.3960 (-1) 


6 


-2.1397 (0) 


-2.0785 (0) 


-1.4856 (0) 


-9.6859 (-1) 


-5.8528 (-1) 


7 


-2.4251 (0) 


-2.4896 (0) 


-1.8690 (0) 


-1.2741 (0) 


-8.2266 (-1) 


8 


-2.6316 (0) 


-2.8676 (0) 


-1.2427 (0) 


-1.5724 (0) 


-1.0515 (0) 


9 


-2.7573 (0) 


-3.2110 (0) 


-2.6069 (0) 


-1.8635 (0) 


-1.2714 (0) 


10 


-2.7981 (0) 


-3.5183 (0) 


-2.9598 (0) 


-2.1471 (0) 


-1.4821 (0) 


11 


-2.7501 (0) 


-3.7882 (0) 


-3.3026 (0) 


-2.4230 (0) 


-1.6834 (0) 


12 


-2.6098 (0) 


-4.0189 (0) 


-3.6343 (0) 


-2.6910 (0) 


-1.8749 (0) 


13 


-2.3736 (0) 


-4.2092 (0) 


-3.9548 (0) 


-2.9510 (0) 


-2.0563 (0) 



" See footnote at end of table. 
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Table 3. ^Coefficients bn(Eo, Z) of the power-series buildup factor formulation (18) derived from the Capo [9] ^ i polynomial coeffi- 
cients according to transformation (d), table 2. — Continued 







bniE,, 


Fe) Z=26 






\^7^:o(Mcv) 












n \^ 


0.5 


0.7 


1.0 


1.5 


2.0 





1.0025 (0) 


9.3280 (-1) 


1.0002 (0) 


1.0162 (0) 


1.0100 (0) 


1 


1.4156 (-1) 


9.3980 (-3) 


2. 1196 (-1) 


3.0483 (-1) 


3.6092 (-1) 


2 


-5.2529 (-1) 


-6.6107 (-1) 


-4.1423 (-1) 


-3.0932 (-1) 


-2.1516 (-1) 


3 


-9.9713 (-1) 


-1.0661 (0) 


-8.7359 (-1) 


-8.2351 (-1) 


-7.1552 (-1) 


4 


-1.2730 (0) 


-1.1932 (0) 


-1.1632 (0) 


-1.2350 (0) 


-1.1374 (0) 


5 


-1.3520 (0) 


-1.0298 (0) 


-1.2726 (0) 


-1.5412 (0) 


-1.4781 (0) 


6 


-1.2332 (0) 


-5.6333 (-1) 


-1.2028 (0) 


-1.7393 (0) 


-1.7348 (0) 


7 


-9.1560 (-1) 


2.1860 (-1) 


-9.4685 (-1) 


-1.8266 (0) 


-1.9048 (0) 


8 


-3.9830 (-1) 


1.3286 (0) 


-5.0014 (-1) 


-1.8004 (0) 


-1.9854 (0) 


9 


3.1965 (-1) 


2.7791 (0) 


1.4217 (-1) 


-1.6581 (0) 


-1.9738 (0) 


10 


1.2392 (0) 


4.5827 (0) 


9.8487 (-1) 


-1.3968 (0) 


-1.8672 (0) 


11 


2.3612 (0) 


6.7518 (0) 


2.0328 (0) 


-1.0139 (0) 


-1.6629 (0) 


12 


3.6868 (0) 


9.2991 (0) 


3.2907 (0) 


-5.0678 (-1) 


-1.3582 (0) 


13 


5.2167 (0) 


1.2237 (1) 


4.7633 (0) 


1.2739 (-1) 


-9.5027 (-1) 


\^Eo(Mev) 












ti \ 


2.5 


3.5 


5.5 


7.5 


9.5 





1.0043 (0) 


9.9882 (-1) 


9.9768 (-1) 


9.9931 (-1) 


1.0011 (0) 


1 


4.1631 (-1) 


5.1544 (-1) 


6.5356 (-1) 


7.3825 (-1) 


7.9417 (-1) 


2 


-1.1129 (-1) 


7.8643 (-2) 


3.4281 (-1) 


5.0375 (-1) 


6.0957 (-1) 


3 


-5.7593 (-1) 


-3.0975 (-1) 


6.5479 (-2) 


2.9451 (-1) 


4.4.503 (-1) 


4 


-9.7500 (-1) 


-6.4792 (-1) 


-1.7841 (-1) 


1.0922 (-1) 


2.9831 (-1) 


5 


-1.3059 (0) 


-9.3405 (-1) 


-3.8883 (-1) 


-5.3419 (-2) 


1.6714 (-1) 


6 


-1.5660 (0) 


-1. 1663 (0) 


-5.6575 (-1) 


-1.9471 (-1) 


4.9286 (-2) 


7 


-1.7526 (0) 


-1.3429 (0) 


-7.0913 (-1) 


-3.1596 (-1) 


-5.7516 (-2) 


8 


-1.8633 (0) 


-1.4620 (0) 


-8.1896 (-1) 


-4.1847 (-1) 


-1.5552 (-1) 


9 


-1.8953 (0) 


-1.5218 (0) 


-8.9519 (-1) 


-5.0353 (-1) 


-2.4696 (-1) 


10 


-1.8461 (0) 


-1.5205 (0) 


-9.3780 (-1) 


-5.7246 (-1) 


-3.3411 (-1) 


11 


-1.7131 (0) 


-1.4562 (0) 


-9.4676 (-1) 


-6.2656 (-1) 


-4.1921 (-1) 


12 


-1.4935 (0) 


-1.3272 (0) 


-9.2203 (-1) 


-6.6712 (-1) 


-5.0452 (-1) 


13 


-1. 1850 (0) 


-1. 1316 (0) 


-8.6360 (-1) 


-6.9545 (-1) 


-5.9228 (-1) 







6n(So, 


Sn) Z=50 






\. £o(Mev) 












71 \ 


0.5 


0.7 


1.0 


1.5 


2.0 





1.0021 (0) 


1.0023 (0) 


1.0014 (0) 


1.0008 (0) 


1.0007 (0) 


1 


4.4293 (-1) 


4.7107 (-1) 


3.8940 (-1) 


3.9964 (-1) 


4.5973 (-1) 


2 


-1.3701 (-1) 


-3.7218 (-2) 


-1.8271 (-1) 


-1.5486 (-1) 


-3.2606 (-2) 


3 


-7.4062 (-1) 


-5.2449 (-1) 


-7.1325 (-1) 


-6.5999 (0) 


-4.7447 (-1) 


4 


-1.3707 (0) 


-9.9267 (-1) 


-1.2005 (0) 


-1.1130 (0) 


-8.6409 (-1) 


5 


-2.0301 (0) 


-1.4437 (0) 


-1.6428 (0) 


-1.5113 (0) 


-1.1997 (0) 


6 


-2.7218 (0) 


-1.8796 (0) 


-2.0384 (0) 


-1.8520 (0) 


-1.4794 (0) 


7 


-3.4484 (0) 


-2.3022 (0) 


-2.3856 (0) 


-2.1326 (0) 


-1.7015 (0) 


8 


-4.2129 (0) 


-2.7134 (0) 


-2.6827 (0) 


-2.3503 (0) 


-1.8642 (0) 


9 


-5.0181 (0) 


-3. 1153 (0) 


-2.9280 (0) 


-2.5023 (0) 


-1.9657 (0) 


10 


-5.8668 (0) 


-3.5097 (0) 


-3.1199 (0) 


-2.5861 (0) 


-2.0042 (0) 


11 


-6.7620 (0) 


-3.8986 (0) 


-3.2566 (0) 


-2.5988 (0) 


-1.9779 (0) 


12 


-7.7063 (0) 


-4.2839 (0) 


-3.3363 (0) 


-2.5378 (0) 


-1.8850 (0) 


13 


-8.7028 (0) 


-4. 6676 (0) 


-3.3576 (0) 


-2.4004 (0) 


-1.7237 (0) 


\ jEJoCMev) 
n N. 












2.5 


3.5 


5.6 


7.5 


9.5 





1.0008 (0) 


1.0010 (0) 


9.9781 (-1) 


9.9514 (-1) 


9.9484 (-1) 


1 


5.1974 (-1) 


6.1402 (-1) 


7.1884 (-1) 


7.8769 (-1) 


8.3584 (-1) 


2 


8.8252 (-2) 


2.7747 (-1) 


4.7050 (-1) 


5.9132 (-1) 


6.7805 (-1) 


3 


-2.9304 (-1) 


-1.0061 (-2) 


2.4198 (-1) 


3.8887 (-1) 


5.0414 (-1) 


4 


-6.2353 (-1) 


-2.4997 (-1) 


1.8862 (-2) 


1.6318 (-1) 


2.9679 (-1) 


5 


-9.0259 (-1) 


-4.4365 (-1) 


-1.9888 (-1) 


-1.0291 (-1) 


3.8674 (-2) 


6 


-1. 1296 (0) 


-5.9248 (-1) 


-4.3285 (-1) 


-4.2656 (-1) 


-2.8753 (-1) 


7 


-1.3039 (0) 


-6.9785 (-1) 


-6.9027 (-1) 


-8.2495 (-1) 


-6.9916 (-1) 


8 


-1.4250 (0) 


-7.6116 (-1) 


-9.8197 (-1) 


-1.3152 (0) 


-1.2135 (0) 


9 


-1.4921 (0) 


-7.8379 (-1) 


-1.3188 (0) 


-1.9146 (0) 


-1.8480 (0) 


10 


-1.5048 (0) 


-7.6713 (-1) 


-1.7204 (0) 


-2.6401 (0) 


-2.6198 (0) 


11 


-1.4622 (0) 


-7.1256 (-1) 


-2. 1708 (0) 


-3.5090 (0) 


-3.5463 (0) 


12 


-1.3639 (0) 


-6.2149 (-1) 


-2.7078 (0) 


-4.5385 (0) 


-4.6449 (0) 


13 


-1.2093 (0) 


-4.9529 (-1) 


-3.3331 (0) 


-5.7457 (0) 


-5.9328 (0) 



*See footnote at end of table. 
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Table 3. "^Coefficients bn(Eo, Z) of the power-series buildup factor formulation {18) derived from the Capo [9] fSi polynomial coeffi^ 
cients according to transformation (d), table 2. — Continued 







6„(^o, W) Z=74 






\^^o(Mev) 












n \ 


0.5 


0.7 


1.0 


1.5 


2.0 





1.0086 (0) 


1.0089 (0) 


1.0032 (0) 


9.9925 (-1) 


9.9840 (-1) 


1 


7.3645 (-1) 


7.1129 (-1) 


5.6840 (-1) 


5.2793 (-1) 


5.6382 (-1) 


2 


4.3004 (-1) 


3.9926 (-1) 


-1.1129 (-1) 


3.9821 (-2) 


1.2485 (-1) 


3 


8.5826 (-2) 


6.7224 (-2) 


-3.7041 (-1) 


-4.6593 (-1) 


-3.2002 (-1) 


4 


-2.9946 (-1) 


-2.9038 (-1) 


-8.7895 (-1) 


-9.9020 (-1) 


-7.7231 (-1) 


5 


-7.2925 (-1) 


-6.7911 (-1) 


-1.4166 (0) 


-1.5338 (0) 


-1.2335 (0) 


6 


-1.2069 (0) 


-1. 1880 (0) 


-1.9855 (0) 


-2.0977 (0) 


-1.7052 (0) 


7 


-1.7358 (0) 


-1.5723 (0) 


-2.5880 (0) 


-2.6827 (0) 


-2. 1888 (0) 


8 


-2.3194 (0) 


-2.0878 (0) 


-3.2264 (0) 


-3.2897 (0) 


-2.6859 (0) 


9 


-2.9610 (0) 


-2.6567 (0) 


-3.9028 (0) 


-3.9195 (0) 


-3.1979 (0) 


10 


-3.6640 (0) 


-3.2846 (0) 


-4.6195 (0) 


-4.5729 (0) 


-3.7265 (0) 


11 


-4.4317 (0) 


-3.9770 (0) 


-5.3788 (0) 


-5.2510 (0) 


-4.2740 (0) 


12 


-5.2677 (0) 


-4.7396 (0) 


-6.1829 (0) 


-5.9545 (0) 


-4.8390 (0) 


13 


-6.1752 (0) 


-5.5767 (0) 


-7.0341 (0) 


-6.6842 (0) 


-5.4261 (0) 


\^o(Mev) 












n \ 


2.5 


3.5 


5.5 


7.5 


9.5 





9.9844 (-1) 


9.9910 (-1) 


9.9375 (-1) 


9.8630 (-1) 


9.8003 (-1) 


1 


6.1013 (-1) 


6.8979 (-1) 


7.6125 (-1) 


7.8912 (-1) 


7.9744 (-1) 


2 


2.2914 (-1) 


4.0578 (-1) 


5.3024 (-1) 


5.6254 (-1) 


5.6596 (-1) 


3 


-1.4706 (-1) 


1.4271 (-1) 


2.8324 (-1) 


2.8009 (-1) 


2.5718 (-1) 


4 


-5.2100 (-1) 


-1.0380 (-1) 


2.7784 (-3) 


-8.4723 (-2) 


-1.5728 (-1) 


5 


-8.9521 (-1) 


-3.3811 (-1) 


-3.2861 (-1) 


-5.5836 (-1) 


-7.0582 (-1) 


6 


-1.2722 (0) 


-5.6459 (-1) 


-7.2839 (-1) 


-1.1673 (0) 


-1.4168 (0) 


7 


-1.6545 (0) 


-7.8762 (-1) 


-1.2140 (0) 


-1.9381 (0) 


-2.3187 (0) 


8 


-2.0446 (0) 


-1.0116 (0) 


-1.8030 (0) 


-2.8971 (0) 


-3.4398 (0) 


9 


-2.4452 (0) 


-1.2408 (0) 


-2.5128 (0) 


-4.0708 (0) 


-4.8086 (0) 


10 


-2.8585 (0) 


-1.4797 (0) 


-3.3608 (0) 


-5.4858 (0) 


-6.4533 (0) 


11 


-3.2873 (0) 


-1.7326 (0) 


-4.3646 (0) 


-7.1685 (0) 


-8.4026 (0) 


12 


-3.7341 (0) 


-2.0039 (0) 


-5.5416 (0) 


-9.1453 (0) 


-1.0685 (1) 


13 


-4.2013 (0) 


-2.2979 (0) 


-6.9093 (0) 


-1.1443 (1) 


-1.3328 (1) 







6n(Eo, Pb) Z=82 






\^Eo(MeT) 












\ 


0.5 


0.7 


1.0 


1.5 


2.0 


n \ 

\ 















9.9993 (-1) 


1.0087 (0) 


1.0159 (0) 


1.0171 (0) 


1.0109 (0) 


1 


7.5580 (-1) 


7.1782 (-1) 


6.7548 (-1) 


6.3718 (-1) 


6.2940 (-1) 


2 


4.7600 (-1) 


3.9601 (-1) 


3.1164 (-1) 


2.4664 (-1) 


2.4906 (-1) 


3 


1.5697 (-1) 


4.0094 (-2) 


-7.5640 (-2) 


-1.5563 (-1) 


-1.3025 (-1) 


4 


-2.0485 (-1) 


-3.5310 (-1) 


-4.9597 (-1) 


-5.7075 (-1) 


-5.0870 (-1) 


5 


-6.1301 (-1) 


-7.8673 (-1) 


-9.4455 (-1) 


-9.9981 (-1) 


-8.8649 (-1) 


6 


-1.0711 (0) 


-1.2639 (0) 


-1.4262 (0) 


-1.4439 (0) 


-1.2638 (0) 


7 


-1.5826 (0) 


-1.7879 (0) 


-1.9433 (0) 


-1.9042 (0) 


-1.6407 (0) 


8 


-2.1512 (0) 


-2.3618 (0) 


-2.4982 (0) 


-2.3818 (0) 


-2.0175 (0) 


9 


-2.7804 (0) 


-2.9887 (0) 


-3.0934 (0) 


-2.8777 (0) 


-2.3944 (0) 


10 


-3.4737 (0) 


-3.6718 (0) 


-3.7314 (0) 


-3.3931 (0) 


-2.7714 (0) 


11 


-4.2347 (0) 


-4.4144 (0) 


-4.4144 (0) 


-3.9291 (0) 


-3. 1488 (0) 


12 


-5.0669 (0) 


-5.2194 (0) 


-5.1449 (0) 


-4.4869 (0) 


-3.5267 (0) 


13 


-5.9709 (0) 


-6.0901 (0) 


-5.9253 (0) 


-5.0674 (0) 


-3.9054 (0) 


\^Eo(Mev) 












n \^ 


2.5 


3.5 


5.5 


7.5 


9.5 





1.0026 (0) 


9.9132 (-1) 


1.0012 (0) 


9.9279 (-1) 


9.6581 (-1) 


1 


6.4235 (-1) 


6.9948 (-1) 


8.0216 (-1) 


8.0445 (-1) 


8.0326 (-1) 


2 


2.9304 (-1) 


4.2933 (-1) 


6.1069 (-1) 


5.8826 (-1) 


5.9791 (-1) 


3 


-4.5410 (-2) 


1.7884 (-1) 


4.1282 (-1) 


3.2014 (-1) 


3.2578 (-1) 


4 


-3.7250 (-1) 


-5.4018 (-2) 


1.9459 (-1) 


-2.4004 (-2) 


-3.7126 (-2) 


5 


-6.8858 (-1) 


-2.7128 (-1) 


-5.7996 (-2) 


-4.6826 (-1) 


-5.1480 (-1) 


6 


-9.9345 (-1) 


-4.5468 (-1) 


-3.5890 (-1) 


-1.0367 (0) 


-1.1312 (0) 


7 


-1.2870 (0) 


-6.6715 (-1) 


-7.2210 (-1) 


-1.7535 (0) 


-1.9104 (0) 


8 


-1.5693 (0) 


-8.4981 (-1) 


-1.1616 (0) 


-2.6426 (0) 


-2.8763 (0) 


9 


-1.8401 (0) 


-1.0250 (0) 


-1.6913 (0) 


-3.7282 (0) 


-4.0529 (0) 


10 


-2.0994 (0) 


-1.1947 (0) 


-2.3252 (0) 


-5.0344 (0) 


-5.4642 (0) 


11 


-2.3472 (0) 


-1.3611 (0) 


-3.0774 (0) 


-6.5852 (0) 


-7.1342 (0) 


12 


-2.5833 (0) 


-1.5260 (0) 


-3.9617 (0) 


-8.4047 (0) 


-9.0869 (0) 


13 


-2.8077 (0) 


-1.6916 (0) 


-4.9921 (0) 


-1.0517 (1) 


-1.1346 (1) 



*See footnote at end of table. 
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Table 3. "^Coefficients b„(Eo, Z) of the power-series buildup factor formulation {18) derived from the Capo [9\ /3i polynomial coeffi- 
dents according to transformation ( d), tab le 2. — Continued 







ftnl-Eo 


U) Z = 92 






\v fi-oCMev) 












n \, 


0.5 


0.7 


1.0 


1.5 


2.0 





1.0107 (0) 


1.0087 (0) 


1.0061 (0) 


1.0033 ( 0) 


1.0017 (0) 


1 


8.3176 (-1) 


7.8986 (-1) 


7.4124 (-1) 


6.9250 (-1) 


6.7550 (-1) 


2 


6.2417 (-1) 


5.4367 (-1) 


4.5187 (-1) 


3.6412 (-1) 


3.3996 (-1) 


3 


3.8503 (-1) 


2.6744 (-1) 


1.3578 (-1) 


1.6510 (-2) 


-6.1238 (-3) 


4 


1.1141 (-1) 


-4.1481 (-2) 


-2.0928 (-1) 


-3.5196 (-1) 


-3.6398 (-1) 


5 


-1.9961 (-1) 


-3.8575 (-1) 


-5.8553 (-1) 


-7.4292 (-1) 


-8.4503 (-1) 


6 


-5.5096 (-1) 


-7.6800 (-1) 


-9.9522 (-1) 


-1.1580 (0) 


-1.1199 (0) 


7 


-9.4557 (-1) 


-1. 1909 (0) 


-1.4406 (0) 


-1.5989 (0) 


-1.5204 (0) 


8 


-1.3864 (0) 


-1.6571 (0) 


-1.9239 (0) 


-2.0671 (0) 


-1.9376 (0) 


9 


-1.8763 (0) 


-2.1692 (0) 


-2.4473 (0) 


-2.5644 (0) 


-2.3726 (0) 


10 


-2.4182 (0) 


-2.7299 (0) 


-3.0131 (0) 


-3.0923 (0) 


-2.8268 (0) 


11 


-3.0151 (0) 


-3.3419 (0) 


-3.6236 (0) 


-3.6526 (0) 


-3.3013 (0) 


12 


-3. 6699 (0) 


-4.0078 (0) 


4.2809 (0) 


-4.2467 (0) 


-3.7974 (0) 


13 


-4.3854 (0) 


-4.7301 (0) 


-4.9874 (0) 


-4.8764 (0) 


-4.3163 (0) 


\ Eo(Mcv) 












\ 


2.5 


3.5 


5.5 


7.5 


9.5 


n \ 












\ 















1.0010 (0) 


1.0006 (0) 


9.9775 (-1) 


9.8835 (-1) 


9.8163 (-1) 


1 


6.8131 (-1) 


7.3122 (-1) 


8.3551 (-1) 


8.2590 (-1) 


8.0789 (-1) 


2 


3. 6080 (-1) 


4.7491 (-1) 


6.8739 (-1) 


6.4374 (-1) 


5.9135 (-1) 


3 


3.8355 (-2) 


2.2977 (-1) 


5.4509 (-1) 


4.2402 (-1) 


3.1040 (-1) 


4 


-2.8712 (-1) 


-6.1556 (-3) 


4.0033 (-1) 


1.4889 (-1) 


-5.6604 (-2) 


5 


-6. 1672 (-1) 


-2.3480 (-1) 


2.4481 (-1) 


-1.9951 (-1) 


-5.3128 (-1) 


6 


-9.5156 (-1) 


-4.5812 (-1) 


7.0266 (-2) 


-6.3902 (-1) 


-1.1353 (0) 


7 


-1.2927 (0) 


-6.7804 (-1) 


-1.3165 (-1) 


-1.1875 (0) 


-1.8902 (0) 


8 


-1.6413 (0) 


-8.9653 (-1) 


-3.6916 (-1) 


-1.8628 (0) 


-2.8176 (0) 


9 


-1.9985 (0) 


-1. 1153 (0) 


-6.5058 (-1) 
-9.8419 (-1) 


-2.6828 (0) 


-3.9393 (0) 


10 


-2.3653 (0) 


-1.3369 (0) 


-3.6653 (0) 


-5.2767 (0) 


11 


-2. 7428 (0) 


-1.5628 (0) 


-1.3783 (0) 


-4.8282 (0) 


-6.8516 (0) 


12 


-3.1321 (0) 


-1.7949 (0) 


-1.8412 (0) 


-6.1893 (0) 


-8.6856 (0) 


13 


-3.5344 (0) 


-2.0353 (0) 


-2.3811 (0) 


-7.7665 (0) 


-1.0800 (1) 



' Tho figures in parentheses in tables 3, 4, and 5 indicate the power of 10 by which the adjacent entry is to be multiplied; 
.,6.4656 (-2) =0.064656. 



6. Utility and Limitations of the Power- 
Series Method 

With respect to use ol' the infinite series (19), 
gamma-ray distributed source problems fall into 
three classes: 

(a) In the region < /jlqxSo .01 , as is often the case 
in an air medium [29], all terms in (19) may be 
neglected except for the zeroth term ((7/47r) go(geoni). 
This term is the complete solution for pure inverse 
square law attenuation as assumed by Meredith 
[29], and depends only on relative rather than abso- 
lute dimensions of the detector-source geometry. 

(b) In the region O.Ol^/io^^l.O the series (19) 
has its greatest usefulness, either with or without 
inclusion of scattered radiation by bn(Eo, Z) coef- 
ficients. This will be demonstrated in the numerical 
examples in section 7. 

(c) In the region 1.0^iulqX<^oo , as is the case in 
fallout shelter calculations, the series (19) is usually 
not useful, and recourse nuist be made to other 
available methods [1, 2]. 

It is interesting to note that the speed of con- 
vergence of (19) is usually enhanced over that for 
(15), rather than the converse. A possible explana- 



tion for this enhancement (with no attempt at a 
proof) is offered as follows: 

In an absorbing and scattering medium the 
point-source buildup factor B:pti{Eo, Z, r) is bounded, 
roughly, by the limits 



l<B,.ri(Eo, Z, r)< expi+fjiov). 



(20) 



The lower bound, Bpti{Eo, Z, r) = l, corresponding to 
jmre absorption (7)^ = in eq (1)), would lead to 
buildup coefficients all equal to unity, i.e., 



bn{Eo, Z) = l, n=0, 1, 2, 



(21) 



in which case (19) would be identical with (15). 
The upper approximate bound, Bfti(Eo,Z, r) 
= exp(+Ato^), corresponding to '^ pure scatter er (pho- 
tons may be scattered, but never absorbed or degrad- 
ed in energy) is the same as simple inverse-square-law 
attenuation, such that the point-source eq (2) would 
reduce to 

D=:{k/4'jr)/r\ (22) 

This buildup-factor behavior can be applied to (19) 
by setting all except the zeroth buildup coeflBcient 
equal to zero, i.e., 
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6,>i(£'o, Z)-0, 



n=l, 2,3, ... . 



(23) 



Hence, in the intermediate region indicated in (20) 
we might expect the observed enhancement. 

7. Numerical Examples 

7.1. Example 1. Rectangular Co-60 Source, Strati- 
fied Barrier of Steel and Water. Comparison 
With Experiment 

There exist little experimental data in the moder- 
ate-penetration region O.OI^mo^^^^I.O, discussed in 
(b) of the preceding section, for simple geometries 
amenable to direct analysis using (19). However, 
in connection with food irradiator design studies [17] 
data have been published for the more complicated 
situation of a two-layer barrier. Here, Donovan 
[30] measured the dose-rate received by a detector 
separated from a rectangular Co-60 source by a 3- 
in. layer of water and varying thicknesses of steel 
cladding as shown in figure 3. These results are 
over the range 0.589 <Ato3:(water+steel)< 1.225, 
hence the}^ provide an excellent check on (19) up to 
the upper limit of its usefulness as discussed in (b) 
of the preceding section. 

The agreement between tlieory (infinite homo- 
geneous medium) and experiment (layered inhomo- 
geneous medium) in the following example must be 
considered somewhat fortuitous, although the posi- 
tion of the detector within the water phantom tends 
to minimize boundary effects. Since interest has 
been expressed in adapting eq (19) for use in nuclear 
engineering curricula, this example is presented in 
great detail, particularly in tlie steps involving- 
dimensional analysis. 

a. Situation 

A 6-in. thick water-equivalent slab phantom is 
separated by an airgap of 2.5 in. from a 50-in. wide 
by 57-in. long rectangular plaque Co-60 source. 
This source, with a uniform activity of 1 curie per 
cm^ of area, is encased in steel cladding of uniform 
thickness. 

b. Problem 

Find the dose-rate (rads [31] in H2O) in the center 
of the phantom as the thickness of steel cladding is 
varied from 0.097 in. to 0.679 in. Use the simplifying 
approximations : 

(1) To take into account the finite thickness of 
the source (Co, Z=27), include half the source 
thickness in the thickness of steel (Fe, Z=26) repre- 
senting '^cladding thickness'^ in the calculation, as 
was done in an alternative analytical study of this 
problem [3]. 

(2) Assume that the Co-60 nuclides emit two 
1.25-Mev photons per disintegration, rather than a 
1.17 Mev and a 1.33 Mev photon as is actuall}^ the 
case. 



CENTERED DETECTOR 



SOLID ANGLE, a 



WATER PHANTOM 




STEEL CLADDING 



Figure 3. Donovan [30] food irradiator rectangular source 
geometry treated in example 1. 



(3) Kalos ([21], page 796) has shown that for a 
composite barrier consisting of a layer of high-Z 
material followed by a layer of low-Z material, the 
resulting buildup factor shows a smooth rise from the 
value for pure high-Z material toward that for pure 
low-Z material, the rate of approach depending on 
the source energy. Indeed, substitution of the linear 
approximation (6) for Bpt,^(Eq, Z, t) in Kalos' 
empirical eq [21], (26.4) gives the linear connecting 
form 

5pTi(£'o, Zi,Z2, ri+rs) 
_ Ato(Zi)ri 



Mo(Zi)ri+/xo(Z2)r2 
, Mo(Z2)r2 



5pTi(£'o,Zi,ri) 

5pTi(fio,22,r2) (24) 



Aio(Z])ri+Mo(Z2)/'2 

in which Zi, n refer to one layer of tlie barrier and 
Z2, r-i refer to the other layer of the barrier. Since 
the buildup for iron and for water are not greatly 
different [32] over the penetration range in this 
problem, the approximate combinatorial form (24) 
will be used. 

c. Plan of Attack 

Evaluate the dose rate (rads in HaO/hr) using (19) 
assuming the medium (a) to be all water: 



7)(H20) = ((7/47r)X:6.(1.25 Mev, H2O) • ^.(geom) • (^o^)^ 

(19a) 
and (b) to be all steel: 



Z?(Fe) = (a/4T)Z;6-.(l-25 Mev, Fe) • g„(geom) • {n^xY 

(19b) 

in which (19a) and (19b) differ only in the values of 
the buildup factor coefficients bniE^, Z). Then use 
(24) to obtain a dose rate between these upper and 
lower bounds by weighting the results from (19a) 
and (19b) by the fractions of mean-free-path distance 
traversed by the primary photons in each material: 

Z^= {(yUoX)H2o/(Mo3^)total}^(H20) 

+ {(MoX)Fe/(Mo^)total}^(Fe). (25) 
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d. Source Strength Constant ((j/47r) 

For an answer from (19) in units of rads/hr per 
curie/cni^ of source strength, tlie factor (o-/47r) is 

(o-/47r) = 1.7005XlO'^XA^(pliotons of energy £'o/disin- 

tegration) X ^o (Mev/photon) X {/x.«(£'o, Z)Ip] (cmVg) 

(26) 

in which the numerical constant 1.7005X10^ is 
obtained from the dimensional factors 

1.7005X10^= (l/47r) (steradians-i)X3.71X10io 

{ (disintegrations/sec)/curie} X1.6X10"^{rads/ 

(Mey/g)}X3.6X10'' (sec/hr) (27) 

For the problem stated above the numerical values 
of the other constants in (26) are 

iV=^ photons/disintegration 

£'o = l»^5 Mev/photon, and 

At.n(1.25 Mev, H2O)/p = 0.0'^96 cmVg.^ 

Substitution of these values in (26) gives a source 
strength constant {gI^tt) of 

((7/47r) = l.*^58XlO^ 

for values of /)(HoO) and D(Fe) from (19a) and (19b) 
in^nits of{rads (H20)/hr} /{curies/cm^} . 

e. Buildup Factor Coefficients h^ (A'o, Z) 

Since reference [9] and hence table 3 of this work 
do not contain data at 1.25 Mev, some interpolation 
is required to obtain the b,X^^o, ^) series coefhcients 
needed in (19a) and (19b). This interj)olation is 
conveniently performed on the '^untraiisformed^' 
l3i{Eo, Z) polynomial coefficients, from which any 
number of bn(Eo, Z) series coefficients inay be gener- 
ated using (d), table 2. Table 4 contains poly- 
nomial coefficients: 

(a) /5,(1.25 Mev, H2O) interpolated from [25],^ 
table 2, and 

(b) I3i(l.25 Mev, Fe) interpolated from [9], table 
7B 



5 Interpolated from [24], table 2. 

6 The coefficients in [25] as in [9] were fitted to tlic calculated results in NYO- 
3075[26]. 



Table 4. Polynomial buildup factor coefficients 13 {{Eq, Z) 
used as input data in examples 1 and 2 



i 


/8,(1.25 Mev, II2O)* 


/3,(1.25Mev, Fe)** 



1 
2 
3 


1.00 (0) 

9.00 (-1) 

1.40 (-1) 

-3.90 (-3) 


1.013 (0) 

7.455 (-1) 

6.055 (-2) 

-5.320 (-4) 



♦Interpolated from reference [25], table 2. 
♦♦Interpolated from reference [9], table 7B. 



Corresponding values of b J 1.25 Mev, H2O) and 
6„(1.25 Mev, Fe), computed using the numbers in 
table 4 as input values for formula (d), tal)le 2, are 
given in table 5 for 0<n<30 and are ready for use 
in (19a) and (19b). 



Table 5. 



Buildup factor and geometry data used in example 
1, series solutions {19a) and {19b). 



The b„(Eo, Z) coefficients, of which 6n(1.25 Mev, H2O) were also used in example 
2, eq (19c), were derived from table 4 using transformation (d), table 2. The 
g-nCgeom) geometry coefficients were obtained from [7], table 3 (and unpublished 
higher coefficients from formulas (57) and (58), as 

5n(Donovan geom) =iqn(a,b) 

using a = (wl2)/h=5, b=(l/2)/h = 5 and the factor of four to apply the detector- 
over-corner data to Donovan's centered-detector geometry. For h=5.5 in., 
these a and b values imply plaque dimensions of 55 by 55 in., rather than 
Donovan's 50 by 57 in. However, the solid angles subtended in each case differ 
by less than 1 percent. 



n 


6„(1.25 Mev, H2O) 


6„(1.25 Mev, Te) 


g„(geom) 





1.000 (0) 


1.013 (0) 


1.090 (1) 


1 


1.000 (-1) 


2.670 (-1) 


-2.953 (1) 


2 


-5.200 (-1) 


-3.574 (-1) 


5.000 (2) 


3 


-8.366 (-D 


-8.575 (-1) 


-6.638 (1) 


4 


-8.264 (-1) 


-1.230 (0) 


7.361 (1) 


5 


-4.660 (-1) 


-1.472 (0) 


-7.022 (1) 


6 


2.680 (-1) 


-1.580 (0) 


5.878 (1) 


7 


1.399 (0) 


-1.551 (0) 


-4.388 (1) 


8 


2.950 (0) 


-1.382 (0) 


2.959 (1) 


9 


4.946 (0) 


-1.069 (0) 


-1.822 (1) 


10 


7.408 (0) 


-6.100 (-1) 


1.032 (1) 


11 


1.036 (I) 


-8.200 (-4) 


-5.424 (0) 


12 


1.383 (1) 


7.613 (-1) 


2. 65S (0) 


13 


1.783 (1) 


1.680 (0) 


-1.220 (0) 


14 


2.240 (1) 


2.757 (0) 


5.274 (-1) 


15 


2.755 (1) 


3.998 (0) 


-2.152 (-1) 


16 


3.330 (1) 


5.404 (0) 


8.322 (-2) 


17 


3.969 (1) 


6.979 (0) 


-3.057 (-2) 


18 


4.673 (1) 


8.726 (0) 


1.069 (-2) 


19 


5.445 (1) 


1.065 (1) 


-3.570 (-3) 


20 


6.288 (1) 


1.275 (1) 


1.140 (-3) 


21 


7.167 (1) 


1.499 (1) 


-3.488 (-4) 


22 


8.192 (1) 


1.750 (1) 


1.024 (-4) 


23 


9.258 (1) 


2.016 (1) 


-2.890 (-5) 


24 


1.040 (2) 


2.300 (1) 


7.852 (-6) 


25 


1.163 (2) 


2.605 (1) 


-2.056 (-G) 


26 


1.294 (2) 


2.929 (1) 


5.197 (-7) 


27 


1.434 (2) 


3.273 (1) 


-1.251 (-7) 


28 


1.583 (2) 


3.577 (1) 


2.998 (-8) 


29 


1.741 (2) 


4.022 (1) 


-6.857 (-9) 


30 


1.908 (2) 


4.429 (1) 


1.520 (-9) 



f. Geometry Coefficients ^n(gcom) 

The parameters describing Donovan^s experi- 
mental geometry, as shown in figure 3 are: 

-?/;= plaque width =50 in. 
Z= plaque length =57 in. 
^z=: distance (height) of detector above plaque center 

= 2.5 in. airgap + 3 in. water 

=5.5 in. 

In the notation of [7] the geometry coefficients 
for this centered-detector situation could be cal- 
culated exactly as 

(/,, (Donovan geometry) =4(/^(a, b) (28) 
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where g^(a, h) are corner-position rectangular-source 
geometry coefficients and 

a=(^/2)/A-:(50/2)/5.5=4.5454 

6=(Z/2)/A=(57/2)/5.5 = 5.18T8 (29) 

However, q^{a, h) values for a=^b, h=b are alreadv 
available in [7] for 0<n<^ and in unpublished 
form for 0<n<35. We can use these tabulated 
coefficients and avoid a two-way interpolation, or an 
exact calculation of (/^(geom) coefficients from 
formulas in [7], by the following argument: 

It has been noted [1-3] that geometry effects scale 
approximately as the solid angle, 12/ subtended by 
the source from the detector for barrier thicknesses 
of the order of one mean-free-path. Hence, since 
the solid angle corresponding to the already tabu- 
lated coefficients differs from the solid angle in 
Donovan's geometry by less than a percent, i.e., 

12(a=5, 5=5)/0(a=4.5454, &-5.1818) 

= 4^o(5, 5)/4:Po(4.5454, S.lsTs) 

= 5.170 steradians/5.132 steradians 

= 1.007 (30) 

the available coefficients for a =5, h = b are~used 
according to (19) and are given as the last column 
of table 5 for 0<n<30. 

g. Penstration Thickness ix^x 

The last factor, {ix^xY, in each series term of (19a) 
and (19b) is for this problem simply the n'i\i power 
of the barrier thickness in mean free paths (abbrevi- 
ated mfp). If this thickness is given in inches, 
the thickness in mfp is 

Mo:r(infp) =^(£^0; Z) (cm^/g) Xp(g/cm') 

X2.54 (cm/in.) Xx(in.) (31) 

where fi(Eo, Z) is the narrow-beam attenuation 
coefficient [23] as used in (2) and p is the density of 
the medium. 

For 1.25 Mev gammas in water we see that 
Mox(mfp)/:r(in.) =0.0634 (cm7g) 

Xl.OO (g/cm3)X2.54 (cm/in.) 
=0.161 mfp/(in. of water) (32) 

and similarly for iron that 

Moa:^(mfp)/x (in.) =0.0533 (cmVg) 

X7.83 (g/cm3)X2.54 (cm/in.) 
=1.06 mfp/ (in. of steel). (33) 



7 This solid angle, ^, is identical with the tabulated quantity po(geom) in [7] 
{i.e., po(a, 6)} and [8] {i.e., :po (p,h)]. 



Thus, for example, 

3 in. water=3 in.X0.161 mfp/in. = 0.483 mfp, 

and 

0.097 in. steel=0.097 in.X1.06 mfp/in. = 0.103 mfp. 

h. Numerical Results, D^ and D 

Donovan's measured values of D for 3 in. of water 
preceded by various thicknesses of steel [30] are 
shown in figure 4, with his estimated ±3 percent 
uncertainty indicated [33]. 

Before evaluating the total dose D according to | 
(19a), (19b), and (25), it is of some interest to first | 
evaluate the unscattered component D^ according 
to (15). We can then note the effect of buildup as 
characterized by the bn(Eo, Z) coefficients. In- 
serting the data from the preceding sections (d), 
(f), and (g) in (15), sums for a variety of /xo^ values 
yield D^ (curve (a) in fig. 4) which is independent 
of the kind of medium and the photon energy 
except for the source strength factor 0-/4 tt. Indeed, 
except for the factor (7/4 tt the data for curve (a) are 
already tabulated in [7], table 4. 

Next we can evaluate /^(HsO) and D(Fe) from 
(19a) and (19b) using the same set of summations 
as for D^ (curve (a), above) except that each series 
term is multiphed by the appropriate buildup 
factor coefficient bniEo,Z) from table 5. These 
sums /^(IlaO) and D(Fe) (for Donovan's geometry 
embedded in infinite media of water and steel, 
respectively) are displayed as curves (b) and (c) 
in figure 4. The circled numbers are a qualitative 
indication of the convergence rate of (19) and are, 
for the indicated calculated points, the number of 
the series term (ri+l, including the zeroth term) 
beyond which the partial sum remains within 0.1 
percent of the sum up to n=30. 

Weighting /^(H^O) and D(Fe) (curves (b) and (c)) 
according to (25) we obtain D, the hybrid curve 
(d), our desired result. 

The calculated values of D in curve (d) can now 
be compared absolutely with Donovan's measured 
values, since the method is free of any arbitrary 
normalization. Agreement is seen to be within 
Donovan's 3 percent estimated uncertainty. 

The results of a more approximate but simpler 
analysis by Moote [3] are shown as curve (e). In 
this analysis, which gives an excellent fit to the data 
in the neighborhood of ijlqX=1, the dose rate D is 
given as the simple product of: 

(a) The source strength constant o-/47r, 

(b) the solid angle ^, indicated in figure 3, sub- 

tended by the source from the detector, 

(c) the exponential integral Ei{ixqx) for the total 

barrier thickness /xoX, and 

(d) the jplane isotropic source buildup factor, 

B^iA{EQ,'H20,iiQx)j for the dose rate at a 
distance /xqX from an infinite plane source 
in water. 
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5 •• EXPERIMENT [30] 



— (a)^ D° Mo-/47r)Z q (geom)-(^ox)^ [7] 

n=0 

(b) D(H^O)= (<T/47r)^Z^bn(Eo,Hj,0) qp(geom)(;iox)" 



(d) D 



30 
h=0 
(MoX)h20D(H20) + (/loX)FeD(Fe) 



(c) D(Fe] = (o-/47r)^2^br^(Eo,Fe)-qn(geom)-(/ioX)" 



(/^o'^)H20 + t/^o'^^Fe 



(e) ? (cr/47r)n(geom)-E,(/xoX)-Bp. ,{Eo,Hp0,/xoX) [3] 




(2 1.25-Mev PHOTONS ASSUMED) 



(n)= NUMBER OF 
^^ TERMS TO SUM 

FOR CONVERGENCE 
TO ^ 0.1 % FROM 
SUM TO n = 30. 



0.5 



I 



= (MoXJH^O 



1.5 



2.0 



+ (MoX)pe 



Figure 4. Results for the Donovan food irradiator rectangular source geometry (example 1), 

Donovan's experimontal points [30, 33] aro comparod with the theoretical results: 

(a) The (los('-ra1(> (-(HiipoiuMit 7)" due to primaries (eq 15). 

(b) The total dosc-ratr /HII2O) assuming an all-water medium (eq 19a). 

(c) The total dose-rate i>(Fe) assuming an all-iron medium (eq 19b). 

(d) A weighted average, D, of curves (b) and (e) according to (25). 

(e) The total dose-rate, D, according to the solid-angle exponential-integral approximation used by Mootc [3]. 



7.2. Example 2. Dose-Rate Profiles Across a Cleared 
Circular Area in a Co-60 Infinite Plane Source 
in Air 

Along the axis of a circular disk source, analytic 
solutions for the flux in terms of exponential integrals 
and the buildup formulations [9, 11-15] are well 
known. However, of -axis solutions including build- 
up data have not, to the author's knowledge, yet ap- 
peared in the literature. The power-series solution 
(19) does apply to this situation, as will be illus- 
trated in this example. 

a. Situation 

A circular area 600 ft in diameter has been cleared 
in an infinite plane isotropic source of 1 curie/ft^ of 
Co-60. The surrounding medium is air (no ground- 
air interface is taken into account). 

b. Problem 

Find dose-rate profiles, relative to the undisturbed 
infinite plane source, at heights 3 ft (0.00626 mfp), 
30 ft (0.0626 mfp) and 300 ft (0.626 mfp) above the 
plane of the source. Buildup data for air are not 



available in parametric form [9, 11-15], but water 
buildup data can be used to good approximation 

[1,2]. 

c. Plan of Attack 

Evaluate the dose-rates (rads/hr in air) at h=?>, 
30, and 300 ft above the infinite plane source [25, 
34] using i(3i(£'o, Z) coeflScients from the polynomial 
buildup formula (10) according to 

Z>PLi(/^)-(<^/47r)27r[Si(Mo/^) 

+ exp(-Mo/iO{/3i+i82+2^3 

+ (^2+2ft)Mo/^+^3(Mo/i')'}] (34) 

in which Ei{iiji) is the ordinary exponential integral 
[35] and values i8^(1.25 Mev, H2O) are already given 
in table 4. Then, using the derived valuesii&;i(1.25 
Mev, H2O) from table 5, evaluate dose-rate profiles 
across a 600-ft diam disk source with the same source 
strength constant (cr/47r) according to (19) 



Z>DisK(p» = (^/47r)i: hn (1.25 Mev, H^O) 

n=0 



q_,{p,h) ■ inJiY (19c) 
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where qn(p,h) geometry coefficients are available in 
[8] for 0<p<10 disk radii and 0.1<A<10 disk radii, 
and in unpublislied form for 0.01<A<0.1. In this 
problem, h assumes the values 0.01, 0.1, and 1.0 disk 
radii while the penetration thickness ^mqIi assumes the 
corresponding values 0.00626, 0.0626, and 0.626 mfp 
given in the problem. 

d. Numerical Results 

Evaluation of (34) using the source strength con- 
stant o-/47r= 12.24 (rads in air/hr)/ (curies/ft^) ^ gives 
infinite-plane dose-rates Dpi,i{h) of 423.7, 247.6, and 
78.83 (rads in air/hr)/(curies/ft^) for the heights 3, 
30, and 300 ft, respectively, as indicated by the hori- 
zontal lines in figure 5. Summation of (19c) for thes3 
three heights and a number of p-values gives the 
dose-rate profiles Ddiqk(p, h) also shown in figure 5. 
The circled numbers, as in example 1, indicate the 
required number of series terms for convergence to 
within 0.1 percent of the stable value. 

Subtraction of Z^disk(p, A) from Z>pli(^) gives the 
desired profiles shown in figure 6, normalized to 100 
percent for each infinite plane source result Z^pli(A), 
as specified in the problem. 



6 This value of (tMtt is equivalent to the Clarke and Buchanan ([34], table] 1; 
value of 14.0 (roentgens/hr)/(curies/ft2) for Co-60. 
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Figure 5. Geometry and intermediate results for example 2. 

The solid, dashed and dotted horizontal lines are infinite-plane-source total dose- 
rates D^-Liih) calculated using (34) for detector heights h=^, 30, and 300 ft 
respectively. The corresponding total dose-rate profiles i)DisK(p,/i) over a 
300-ft radius disk source were calculated using (19c). 
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Figure 6. Relative dose-rate profiles at heights 3, 30, and 300 ft 
above a 300-ft radius cleared circular area in an infinite plane 
Co-60 source in an infinite air medium {example 2). 

These profiles were obtained by subtracting the disk-source profiles Z>disk(p,/i) 
in figure 5 from the corresponding infinite plane source dose-rates Dvu.{h), 
dividing by DpLi(h) and multiplying by 100 percent. 

Caution must be taken in applying these results to 
the study of radiation fields from fallout deposition 
on the ground, since the ground-air interface has 
been shown by theoretical [36] and experimental [37, 
38] studies to affect the scattered component by as 
much as a factor of three. Qualitatively, this would 
have the greatest effect on the h=3 ft profile in 
figure 5, ''squaring the shoulders'' of the curve and 
depressing the tails on either side. In figure 6 the 
3-ft profile would more closely resemble a '^square 
weff,'' but the 30-ft and 300-ft profiles would not be 
appreciably changed. 

8. Discussion 

The power-series formulation for buildup data is 
seen to have two main limitations: 

(a) Convergence of the primary-flux series (15), 
and hence the total-fiux series (19), is poor for 
combinations of thick barriers (ijlqX^I) and detector- 
source separations much less than the source dimen- 
sions (i.e., a, 6>>1 or h<C<Cl). In such cases the 
geometry closely resembles the infinite plane source 
geometry for which the series (15) for D^ becomes, 
in the limit, the series representation ((5), [35]) for 
the exponential integral Ei(iiox). The latter series 
converges absolutely for 0<^^to^<^oo^ but for 1< 
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Mo^<C °° the convergence rate is rather slow for prac- 
tical use. In these situations it is more economical 
to revert to numerical integration, or to use an 
infinite plane source analytical formulation, such as 
(34) in example 2, as a close approximation. 

(b) Boundary ell'ects such as tlie density-interface 
effect [36-38] are neglected, but this limitation is 
inherent in any buildup factor formulation. Prob- 
lems in which such effects are important are best 
treated by Monte Carlo techniques at present [36]. 

On the other hand, the series formulation provides 
some advantages : 

(a) The three-way separation variables [(1) unique 
properties of the medium — bn{Eo,Z), (2) geometry— 
qn (geom) and (3) barrier thickness — -ixqx] gives both 
economy of data tabulation and flexibility of ap- 
plication. A given set of buildup data can be a])phed 
to a variety of distributed source geometries, or 
vice versa, and at the same time the barrier thickness 
/jlqX can be varied in a trivial fashion. 

(b) In some simple idealized geometries the series 
(19) can provide exact answers, especially in the 
region of small barrier thickness (aio'^<1). These 
answers, in addition to available experimental data, 
can then be used to check the validity of simpler 
engineering-type approximate methods (e.g., [3J) in 
this region. 



The autlior tlianks R. L. Bach for computing 
table 3, R. J. Herbold for com])uting the additional 
qnigeom.) coefficieiits needed for example 2, and 
M. J. Berger for his interest and helpful criticisms 
regarding this effort. 
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